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Abstract: The fact that nuclei have diffuse surfaces (rather than being simple spheres) has 
dramatic consequences on the interpretation of the RHIC heavy-ion data. The effect is quite 
small (but not negligible) for central collisions, but gets increasingly important with decreasing 
centrality. One may actually divide the collision zone into a central part ("core"), with expected 
high energy densities, and a peripheral part ("corona"), with smaller energy densities, more like in 
pp or pA collisions. We will discuss that many complicated "features" observed at RHIC become 
almost trivial after subtracting the corona background. We are focussing on AuAu collisions at 
200 GeV. 



Everybody having worked on pA scattering 
knows about the importance of the nuclear dif- 
fuseness. Even though the proton suffers on 
the average up to 6 collisions when it traverses 
a big nucleus, the most likely situation is still 
just one interaction, due to the surface diffuse- 
ness. 

In nuclear collisions, the surface effect is as 
well present, and also very important. The pe- 
ripheral nucleons of either nucleus essentially 
perform independent pp or pA-like interac- 
tions, with a very different particle production 
compared to the high density central part. For 
certain observables, this "background" contri- 
bution completely spoils the "signal", and to 
properly interpret RHIC data, we need to sub- 
tract this background. 

In order to get quantitative results, we need 
a simulation tool, and here we take EPOS fl], 
which has proven to work very well for pp and 
dAu collisions at RHIC. It is very important 
to understand that the main results of this 
paper do not depend on whether or not the 
model treats the high density part 100% 
correctly. The crucial point is that the 
model describes pp and pAu to a high pre- 
cision, so we can safely subtract the periph- 
eral low density part. 



EPOS is a parton model, so in case of a 
AuAu collision there are many binary interac- 
tions, creating partons, which then hadronize 
employing a phenomenological procedure call 
string fragmentation. Here, we modify the pro- 
cedure: we have a look at the situation at an 
early proper time tq, long before the hadrons 
are formed: we distinguish between string seg- 
ments in dense areas (more than pa segments 
per unit proper volume), from those in low den- 
sity areas. In the following, we will use ro = 1 
and po = 1 fm^^. We refer to high density areas 
as core, and to low density areas as corona. 

In figure. ^ we show two examples (ran- 
domly chosen) of semi-peripheral (40-50%) 
AuAu collisions at 200 GeV (cms), simulated 
with EPOS. We observe large fluctuations event 
by event, simply reflecting the randomness of 
the positions of binary nucleon-nucleon colli- 
sions. There is an important contribution from 
the low density area, contributing roughly 20% 
to the final particle production. But much more 
importantly, as discussed later, the importance 
of this contribution depends strongly on parti- 
cle type and transverse momentum. For cen- 
tral collisions, the low density contribution is 
obviously less important, for more peripheral 
collisions this contribution will even dominate. 

How do these low density contributions inter - 
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Figure 1 : Two Monte Carlo realizations of semi- 
peripheral (40-50%) AuAu collisions at 200 GeV 
(cms) . We show in red string segments in high 
density areas (core), and in green the string 
segments in low density environments (corona). 
The circles are put in just to guide the eye: they 
represent the two nuclei in hard sphere approx- 
imation. We consider a projection of segments 
within z = ±0.8 fm to the transverse plane (x,y). 



act with the expanding core? Well, even a sys- 
tem of noninteracting particles expands, with 
the velocity of light (reflecting the outward mov- 
ing particles ). Inward moving particles may 
be absorbed by the core, on the other hand 
the core edges start to hadronize at the same 
time, with a good chance that early hadroniza- 
tion and absorption compensate each other. So 
we ignore any interaction for the moment. But 
even if part of the low density contribution will 
be absorbed, there will be a sizable effect. 

In order to make a quantitative statement, 
we adopt the following strategy: the low den- 
sity part will be treated using the usual EPOS 
particle production which has proven to be very 
successful in pp and dAu scattering (the pe- 
ripheral interactions are essentially pp or pA 
scatterings). For the high density part, we sim- 
ply try to parameterize particle production, in 
the most simple way possible (It is not at all 
our aim to provide a microscopic description 
of this part). Suppose we find such a sim- 
ple parameterization of the core contribution, 
such that the total contribution reproduces all 
the relevant low and intermediate -pt data, then 
our core parameterization represents in fact the 
data after "background subtraction", and that 
is what we are really interested in! 

In practice, we first divide the EPOS string 
segments into core and corona contribution, as 
discussed earlier (apart of the density, there is 
another condition: only segments with trans- 
verse momenta less than 4 GeV contribute 
to the core, the others escape freely, no Jet 
quenching). We then consider the core con- 
tributions more closely, in longitudinal slices, 
characterized by some range in t] ~ 0.51n(t + 
z)l(t — z). Since string segments show a 
Bjorken-fluid-like behavior, the particles in a 
segments around rj move with rapidities close 
to T]. Connected core regions in a given segment 
are considered to be clusters, whose energy 
and fiavor content are complete determined by 
the corresponding string segments. Clusters 
are then considered to be collectively expand- 
ing: Bjorken-like in longitudinal direction, with 
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in addition some transverse expansion. We 
assume that the clusters hadronize at some 
given energy density Ehadi . having acquired at 
that moment a collective radial flow, with a lin- 
ear radial rapidity profile from inside to out- 
side, characterized by the maximal radial ra- 
pidity j/rad- In addition, we impose an azimuthal 
asjonmetry, by multiplying the x and y compo- 
nent of the flow four -velocity with 1 + e f^cc and 
1 — e /ccc. where e is the the initial spacial eccen- 
tricity, e = {{y'^ - x^)/y'^ +x^)), and /ccc a pa- 
rameter. By imposing radial flow, we have to 
rescale the cluster mass as 

M M X 0.5 y^ad/ (yrad sinhj/rad - coshyrad + 1), 

in order to conserve energy. Hadroniza- 
tion then occurs according to covariant phase 
space, which means that the probability dP of 
a given final state of n hadrons is given as 

1 " ^3 
species a i—1 

with Pi — {Ei,pl) being the four -momentum of 
the i-th hadron, gi its degeneracy, and /; its 
quark flavor content [u — u,d — d...]. There is a 
factor Si = 7s for each strange particle (sign 
plus for a baiyon, sign minus for a meson), 
with 7s being a parameter. The number 
counts the number of hadrons of species a. E 
is the total energy of the cluster in its cms, W 
is the cluster proper volume. The whole pro- 
cedure perfectly conserves energy, momentum, 
and flavors (microcanonical procedure). 

So the core definition and its hadronization 
are parameterized in terms of 6 global parame- 
ters: 



To 


1 fm 


core formation time 


Po 


1 fm-'* 


core formation density 


£hadr 


0.22 ^ 


hadr. energy density 


yrad 


0.83 


max. radial flow rapidity 


/ccc 


0.5 


eccentricity coefflcient 


Is 


1.3 


hadronization factor 



The flnal results are insensitive to variations 
of To, even changes as big as a factor of 2 do 



not affect the results at all. This is a nice fea- 
ture, indicating that the very details of the ini- 
tial state do not matter so much. We call these 
parameters "global", since they account for all 
observables at all possible different centralities 
at RHIC. In the following, we are going to dis- 
cuss results, all obtained with the above set of 
parameters. 

All the discussion of RHIC data will be based 
on the interplay between core and corona con- 
tributions. To get some feeling, we flrst com- 
pare in flg. 121 the core contribution correspond- 
ing to a central (0-5%) AuAu collision (which 
means purely statistical hadronization, with 
flow) with pp scattering (which is qualitatively 
very similar to the corona contribution). We 
plot nit spectra of pions, kaons, protons, and 
lambdas, the nuclear spectra are divided by 
the number of binary collisions (according to 
Glauber). We observe several remarkable fea- 
tures: the shapes of the different pp spectra 
are not so different among each each other, 
there is much more species dependence in the 
core spectra, since the heavier particles acquire 
large transverse momenta due to the flow ef- 
fect. The second main observation concerns 
the yields, in particular at intermediate val- 
ues of mt — m: the yields for the different pp 
contributions are much wider spread than the 
core contributions; in particular, pion produc- 
tion is suppressed in the core hadronization 
compared to pp, whereas lambda production 
is favored. All this is quite trivial, but several 
"mysteries" discussed in the literature (and to 
be discussed later in this paper) are just due to 
this. 

Let us now compare core and corona contri- 
butions for different centralities in AuAu col- 
fisions at 200 GeV. In flg. |3l we plot the rel- 
ative contribution of the core (relative to the 
complete spectrum, core + corona) as a func- 
tion of mt~m, for different particle species. For 
central collisions, the core contribution domi- 
nates largely (around 90%), whereas for semi- 
central collisions (40-50%) and even more for 
peripheral collisions the core contribution de- 



3 




2.5 3 
mt-m 

Figure 2: Invariant yields l/2Trmtdn/dydmt of 
pions (red), kaons (blue), protons (green), and 
lambdas (yellow), for the core contribution cor- 
responding to a central (0-5%) AuAu colli- 
sion (full lines) and proton-proton scattering 
(dashed). 

creases, giving more and more space for the 
corona part. Apart of these general statements, 
the precise mt dependence of the relative weight 
of core versus corona depends on the particle 
type, and can be easily understood by inspect- 
ing figure 121 since the corona contribution is up 
to a factor very close to pp. 

We are now ready to investigate RHIC data. 
In fig. 12 we plot the centrality dependence of 
the particle yield per participant (per unit of ra- 
pidity), for 7r+, K+, and p, the data together 
with the full calculation, but also indicating the 
core contribution. The complete calculation fol- 
lows quite closely the data. Whereas central 
collisions are clearly core dominated, the core 
contributes less and less with decreasing cen- 
trality. Similar results are obtained for tt^, , 
and p, and also lambdas and xis. 

Next we consider particle ratios, as a func- 
tion of centrality. In fig. [51 we show the ra- 
tios of different particles, with respect to pi- 
ons. Whereas the complete contribution (as the 
data) show a strong centrality dependence, the 
rations are practically flat for the core contri- 
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Figure 3: The relative contribution of the core 
(core/(core+corona)) as a function of the trans- 
verse mass for different centralities (0-5%: red, 
40-50%: blue, 70-80%: green). Upper figure: 
pions (full) and kaons (dashed). Lower figure: 
protons (full) and lambdas (dashed). 
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Figure 4: Rapidity density dn/dy per partic- 
ipant as a function of the number of partici- 
pants, for 7r+ (red). A' + (blue), and p (green). We 
show data (points) t2J together with the full cal- 
culation (core + corona, fuU line) and just the 
core part (dashed) . 



< 
< 



3.5 F 




Figure 6: Nuclear modiflcation factors in cen- 
tral AuAu collisions at 200 GeV. Lines are full 
calculations, symbols represent data fT,^. We 
show results for pions (red; circles), kaons 
(blue; squares), protons (green; triangles), and 
lambdas (yellow; inverted triangles). 
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Figure 5: Particle ratios as a function of cen- 
trality: /i+/7r+(red), A'^/tt^ (green), p/7r+ (blue), 
p/7r~(cyan), A/7r^(gray), S^/tt^. Complete cal- 
culation (full) and just ratio of the core contri- 
butions (dotted). 

buttons, apart of some decrease for very small 
participant numbers. 
So our first important conclusion: af- 



ter subtracting the "corona background", 
the interesting part, the core contribution, 
shows an extremely simple behavior: there 
is no centrality dependence, the systems are 
simply changing in size (and the participant 
number is certainly not a good measure of 
the volume of the core part, this is why 
the overall multiplicities per participant de- 
crease with decreasing centrality). 

Lets us come to pt or nit spectra. We checked 
all available low and intermediate pt data (pi- 
ons, kaons, protons, lambdas, xis), and our 
combined approach (core + corona) describes 
well the data (better than the differences be- 
tween STAR and PHENIX results). Lacking 
space, we just discuss a (typical) example: the 
nuclear modification factor (AA/pp/number of 
collisions), for pions, kaons, protons, and 
lambdas in central AuAu collisions at 200 GeV, 
see fig. |6l For understanding these curves, we 
simply have a look at fig. |2l where we com- 
pare the core contributions from AuAu (divided 
by the number of binary collisions) with pp. 
Since for very central collisions the core domi- 
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nates largely, the ratio of core to pp (the solid 
line divided by the dotted one, in fig. |2j corre- 
sponds to the nuclear modification factor. We 
discussed already earlier the very different be- 
havior of the core spectra (phase space decay) 
compared to the pp spectra (string decay): pi- 
ons are suppressed, whereas heavier particles 
are favored. Or better to say it the other way 
round: the production of baryons compared to 
mesons is much more suppressed in string de- 
cays than in statistical hadronization. 

So what we observe here, is nothing but 
the very different behavior of statistical 
hadronization (plus flow) on one hand, and 
string fragmentation on the other hand. 
This completely statistical behavior indi- 
cates that the low pt partons do not suffer 
energy loss, they get completely absorbed in 
the core matter. 

The Rep modification factors (central over 
peripheral) are much less extreme than RaA' 
since peripheral AuAu collisions are a mix- 
ture of core and corona (the latter one being 
pp-like), so a big part of the effect seen in 
Raa is simply washed out. 

Let us finally discuss elliptical flow, shown 
in fig. [71 We understand the results in the 
following way: the pion curve seems to satu- 
rate at high Pt , which is here simply due to the 
fact that with increasing pt the continuously 
increasing core curve is more and more "con- 
taminated" by corona contributions. For the 
lambdas, the effect is much smaller, since the 
corona contributions are smaller, as seen from 
fig. |3l Eventually, the lambda curve will also 
saturate, but at larger pt. 

To summarize: we have discussed the infiu- 
ence of the corona contribution (occurring in 
the periphery of nuclear collisions) in AuAu col- 
lisions at RHIC. Our analysis is based on a 
model which works excellently for pp and pA, 
together with a very simple parameterization of 
the central (core) part. The fact that this sim- 
ple treatment works, indicates that the part we 
are really interested in, the core, shows a very 
simple behavior. For example, contrary to the 




0.5 1 1.5 2 2.5 3 3.5 4 

pt 



Figure 7: Elliptical fiow for pions (red) and 
lambdas (green), as a function of transverse 
momentum. The points are data jSj |3, lines 
are calculations. The dotted lines represent 
only the core contribution, the full lines are the 
complete contributions, core + corona. 

general believe, there seems to be no central- 
ity dependence of particle production. Just the 
volume changes. 

We do not make any attempt here to ex- 
plain these very interesting data, the only pur- 
pose here is to separate the interesting part 
(core) from the contamination (corona) . We also 
did not make any efforts to optimize the fits, 
actually most parameters are essentially first 
guesses. To get more precision one need to en- 
ter into a more technical discussion about for 
example the feed-down correction procedures 
in the different experiments. 
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